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1.  INTRODUCTION 


The  purpose  of  this  report  is  to  describe  the  thennosonde  instrument  used  to  obtain  the 
vertical  profile  of  the  refractive  index  structure  constant,  C„2.  The  balloon-borne 
thennosonde  instrument  is  a  package  originally  developed  by  NASA,  then  improved  over 
the  years  by  the  Air  Force  Geophysics  Laboratory,  now  AFRL/VSBL.  The  instrument  has 
evolved  through  the  years,  but  basically  it  involves  measurements  of  the  temperature 
difference  across  a  fixed  horizontal  distance,  typically  lm.  This  is  done  while  the  balloon 
ascends  at  an  approximately  constant  vertical  velocity,  5  m/s.  This  temperature  difference 
is  used  to  obtain  the  temperature  structure  constant,  Ct  .  The  refractive  index  structure 
constant,  Cn,  is  obtained  from  Cf  by  using  of  the  Dale-Gladstone  equation. 

In  the  years  since  the  last  thermosonde  technical  report1  in  1982,  calibration  procedures 
have  gone  through  several  changes.  The  evolution  began  with  a  change  in  radiosonde 
telemetry  method  from  analog  frequency  changes  to  digital  “counts”.  Lower  instrument 
noise  has  allowed  the  dual  gain  calibration  curves  to  be  replaced  by  a  single  gain 
calibration  curve.  The  multiple  segment  calibration  relationship  was  replaced  by  a  single 
straight  line.  A  new,  lower  signal  calibration  input  point  has  been  added  to  the  calibration 
source.  As  instrument  noise  has  reduced  and  component  and  assembly  quality  has 
improved,  individual  instrument  calibration  curves  have  been  replaced  by  a  single 
universal  calibration  curve.  The  approach  to  handling  instrumentation  noise  has  also 
changed. 

The  purpose  of  this  report  is  to  document  the  current  thennosonde  configuration.  An 
important  section  is  a  description  of  the  current  calibration  methodology  and  the  new 
universal  calibration  curve.  Also  documented  are  the  uncertainties  associated  with  the 
current  hardware  and  new  calibration  procedures. 


2.  CN2  MEASUREMENT  PRINCIPLES 

The  basis  for  the  thennosonde  measurement  is  the  analysis  of  random  processes.  To  the 
extent  that  the  fluctuations  of  a  variable,  Y,  in  the  atmosphere  can  be  considered  a 
homogeneous  and  isotropic  random  process,  one  can  pass  two  sensors  through  the 
atmosphere  a  distance,  r,  apart  and  generate  the  structure  function,  Dy  (r  )  where2 

A(r)  =  ([T(x)-T(x  +  r)]2).  (1) 
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A.  N.  Kolmogorov3  postulated  that  for  velocity  fluctuations  in  a  fully  developed  turbulent 
flow,  the  velocity  structure  function,  Dv(r),  would  have  the  following  equivalent  relation 
to  the  distance  in  the  structure  function: 

Dv{r)  =  Cv2r2'3  (2) 

where  the  constant  C2  is  known  as  the  velocity  structure  constant.  If  one  were  to 
observe  the  velocity  power  spectral  density,  Sn.(K ) ,  it  would  have  the  equivalent  relation: 

Srr  ( K )  =  0.25 C2vK~5/3  (3) 

where  K  is  the  wave  number  of  the  eddies. 

The  random  turbulence  field  was  assumed  to  be  locally  homogeneous  and  isotropic  for 
scales  less  than  the  largest  eddies,  which  are  the  turbulent  kinetic  energy  source  (the  outer 
scale,  L0),  but  larger  than  the  smallest  eddies,  which  are  the  sizes  of  the  viscous  energy 
dissipation  (the  inner  scale,  lo).  The  range  between  the  scales,  where  the  turbulent  kinetic 
energy  is  passed  down  from  the  larger  eddies  to  the  smaller  by  inertial  processes,  is 
known  as  the  inertial  range.  While  the  condition  of  inertial  range  “equilibrium 
turbulence”  defined  above  is  not  always  present  in  the  atmosphere,  functional 
relationships  with  a  power  near  2/3  are  surprisingly  prevalent. 

Corrsin  ,  Yaglom,  and  Obukhov  ’  independently  extended  the  Kolmogorov  theory  to 
passive  scalars,  and  showed  that  the  structure  function  of  conservative  passive  scalars, 
such  as  temperature  and  refractive  index,  should  also  follow  the  2/3  power  law  (Eq.  2). 

In  this  report,  it  is  the  refractive  index  structure  constant,  C2 ,  that  is  of  interest.  It  is 
obtained  from  the  temperature  structure  constant,  C] ,  via  the  Dale-Gladstone  relation,  i.e., 


( 

79xl0"6 

v 


In  the  above  expression,  P  is  the  atmospheric  pressure  in  millibars  and  T  is  the 
temperature  in  degrees  Kelvin. 


(5) 


3.  DESCRIPTION  OF  THE  INSTRUMENT 

The  fine  resistance  wires  used  in  the  thermosonde  have  a  time  constant  less  than  2 
ms  at  sea  level.  This  pennits  the  measurement  of  rapidly  varying  temperature 
fluctuations.  Low  constant  current  passes  through  the  wires,  which  operate  as  low- 
overheat  resistance  wire  detectors.  The  two  resistance  probes  are  legs  of  a  Wheatstone 
bridge,  which  generates  a  voltage  proportional  to  the  resistance  difference.  The 
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instrumentation  must  be  sensitive  enough  to  detect  temperature  differences  in  milli- 
degrees,  Kelvin.  The  fluctuating  difference  voltage  is  amplified,  filtered,  and  passed  to 
an  onboard  analog  processor  that  computes  the  root  mean  square  (RMS)  of  the 
fluctuating  temperature  difference.  The  time  constant  of  the  RMS  chip  is  set  to  3.7s.  The 
thermosonde  instrument  is  attached  to  a  standard  meteorological  radiosonde  measuring 
the  ambient  temperature,  pressure,  humidity  (Figure  1).  The  RMS  chip  output  is  sampled 
every  1  to  2  seconds  (depending  on  the  type  of  radiosonde)  by  a  spare  channel  of  the 
radiosonde  and  relayed  to  the  ground  with  the  meteorological  data.  The  radiosonde  also 
uses  either  GPS  or  the  Loran  navigation  system,  from  which  position  and  wind  velocity 
are  determined  and  relayed  down  to  the  ground  station.  The  ground  station  has  special 
software  to  extract  the  spare  channel  data.  A  functional  schematic  of  the  thermosonde 
system  is  shown  in  Figure  2. 

In  order  to  understand  the  response  of  the  thermosonde,  an  engineering  model  was 
assembled  and  analyzed  in  a  computer  simulation.  The  basis  of  the  engineering  model 
are  two  AFGL  technical  reports1'6,  the  NASA  technical  note  by  the  thermosonde 
developer,  Jack  Bufton  ,  and  recent  inhouse  engineering  work  .  The  simulation  model 
uses  the  software  package,  Advanced  Continuous  Simulation  Language  (ACSL)9.  A 
graphical  model  of  the  system  is  shown  in  Figure  3. 


Figure  1 .  Two  views  of  the  thermosonde  showing  probes,  beam  and  radiosonde. 
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Figure  2.  Functional  Schematic  of  the  thermosonde  integrated  with  the  VIZ  radiosonde 


Figure  2.  A  graphical  representation  of  the  thennosonde  electronics  in  ACSL  simulation. 
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Referring  to  the  graphical  representation  in  Figure  3,  there  are  some  other  technical  details 
that  give  a  more  complete  description  of  the  instrument.  The  calibration  box  is  attached  to 
the  bridge,  in  place  of  the  probes,  using  a  procedure  described  below. 

Since  there  can  be  some  drift  of  resistance  difference  as  the  instrument  changes 
temperature  as  it  is  ascending  through  the  atmosphere,  there  is  a  bridge  rebalancing  circuit 
with  a  time  constant  of  20s. 

The  overall  gain  of  the  instrument  is  approximately  600,000,  which  is  possible  due  to  an 
aspect  of  the  circuit  that  is  not  part  of  the  simulation,  synchronous  detection.  The  bridge 
and  the  initial  part  of  the  circuit,  is  actually  powered  by  alternating  current  at  3kHz.  The 
desired  signal  is  amplitude  modulated  on  this  3kHz  signal.  After  most  of  the 
amplification,  a  synchronous  detector  strips  out  the  correct  frequency. 

After  low  and  high  pass  filtering,  the  signal  is  passed  to  a  root-mean-square  integrated 
circuit  or  “RMS  Chip”.  This  does  not  perform  an  exact  mathematical  equivalent  of  the 
RMS  operation.  Whereas  the  chip  is  often  represented  by  a  first  order  integrator  with  a 
fixed  time  constant,  the  circuit  actually  gives  different  time  constants  for  the  rise  and  fall 
of  the  input  signal  .  The  output  of  the  RMS  chip  is  sampled  and  passed  to  the  radiosonde 
for  transmission  to  the  ground  station. 


4.  CALIBRATION 

This  section  documents  the  current  calibration  process  used  for  thermosondes 


4.1  Telemetry  Calibration 

The  initial  calibration  transfonnation  is  from  received  digital  counts  to  RMS  voltage. 
Two  brands  of  radiosondes  have  been  used  through  the  period,  VIZ  (now  Lockheed 
Martin  -  Sippican)  and  Vaisala.  The  VIZ  sondes  were  individually  calibrated  prior  to 
flight,  and  have  been  found  to  have  very  little  variation  from  sonde  to  sonde. 

The  Vaisala  sondes  have  individual  calibration  constants  for  the  interface  box  through 
which  the  thermosonde  electronics  are  connected  to  the  radiosonde.  There  can  be  large 
variations  from  box  to  box.  The  calibrations  are  checked  periodically  and  found  to  be 
accurate.  The  calibration  constants  are  part  of  the  input  into  the  data  reduction  program. 
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4.2  Instrument  Noise 


When  calibrating  the  thennosonde  instrument,  one  setting  is  an  open  circuit,  which 
equates  to  a  zero  input  signal.  With  this  input,  there  is  an  output  voltage  from  the  system, 
which  is  interpreted  as  the  instrument  noise  floor.  In  the  laboratory,  this  noise  floor  is 
0.019  to  0.025  Volts.  When  the  zero-signal  input  to  the  RMS  chip  was  sampled  in  the 
laboratory,  it  was  found  to  be  a  Gaussian  distribution  with  a  mean  of  0  and  a  standard 
deviation  of  approximately  the  output  noise  floor  of  the  chip.  After  the  instrument  is 
released,  and  ascending  through  the  atmosphere  behind  a  balloon,  the  minimum  voltage  is 
usually  a  few  mV  lower  than  the  noise  floor  observed  in  the  laboratory.  This  is  not 
considered  unusual  since  electronic  noise  is  higher  inside  buildings  due  to  light  bulbs  and 
other  electrical  appliances  and  instruments.  The  apparent  flight  noise  floor  is  often  not 
constant,  but  may  slowly  drift  lower  or  higher  as  the  flight  progresses.  A  representative 
value  of  the  inflight  noise  is  recorded  and  entered  into  the  data  reduction  program. 

The  current  method  of  handling  instrument  noise  is  to  assume  that  the  noise  is 
independent  of  input  signal  and  uncorrelated  with  that  signal,  and  that  the  noise  amplitude 
is  only  slowly  varying  with  environmental  conditions.  The  signal  into  the  RMS  chip  is 
then: 


K  =  K  +  K  (6) 

Where  the  subscripts  represent  input,  signal,  and  noise  respectively.  The  output  of  the 
RMS  chip,  Vrms ,  can  be  represented  by  the  following  equation: 


V  = 

rms 


k(K+K)2J‘T 


(7) 


where  r  is  the  time  constant  of  the  RMS  chip  which  is  long  compared  to  the  fluctuations 
of  the  turbulent  signal  and  the  noise.  Performing  the  square,  and  recognizing  that  the 
integral  is  the  time  average  of  the  quantity  over  r  ,  we  obtain  the  result: 

vL={  K2)+2{KK)+{K2}  (8) 

where  the  brackets  indicate  the  time  average.  Before  entering  the  RMS  chip,  the  noise  is 
equally  distributed  in  the  positive  and  negative  direction.  In  addition,  the  bridge 
balancing  circuit  insures  that  the  input  signal  is  also  balanced  about  a  mean  of  zero  volts, 
to  within  the  time  constant  of  the  balancing  circuit  (~20s).  Further,  we  assume  that  there 
is  no  correlation  between  signal  and  noise.  The  result  is  that  the  middle,  cross  product 
tenn  averages  out  to  nearly  zero,  leaving  us  with  the  equation: 

VL={K2)  +  {V,2)  (9) 


6 


When  the  input  signal  is  zero,  there  is  a  small  output  from  the  chip.  This  is  interpreted  as 
the  noise.  The  RMS  chip  has  converted  the  oscillating  circuit  noise  to  a  positive  value, 
which  is  averaged  over  time,  and,  in  fact,  does  not  vary  significantly  over  a  reasonable 
period  of  observation  on  the  bench.  Since  we  are  interested  in  the  amplified  signal  from 
the  bridge,  Vs ,  the  RMS  chip  output  voltage  is  always  corrected  for  the  noise.  In 

calibration,  the  noise  is  determined  when  no  resistance  change  is  introduced  from  the 
calibration  box,  and  each  subsequent  calibration  output  is  corrected  using  the  formula: 

ao> 

Where  we  have  dropped  the  time  average  brackets,  and  each  voltage  is  assumed  to  be  a 
time  averaged  value.  This  procedure  is  consistent  with  guidance  in  correcting  for  noise 
in  AC  voltmeters1  \  This  correction  is  insignificant  for  higher  values  of  signal,  but 
becomes  significant  as  the  signal  strength  decays  to  low  values. 

4.3  Calibration  Theory 

It  is  assumed  that  the  change  in  resistance  of  a  material  with  a  change  in  temperature  is 
nearly  linear,  and  for  modest  temperature  differences,  it  can  be  expressed  as: 

R(T)  =  R0  [1 + a  ( r  -  r0 )]  (ii) 

where  the  subscript,  0,  denotes  a  reference  temperature,  usually  chosen  to  be  0°C.  Where 
a  is  the  temperature  coefficient  of  resistivity,  and  it  has  a  value  of  0.00375K’1  for  the 
thennosonde  probe  wire. 


In  the  thennosonde,  the  Wheatstone  bridge  detects  a  difference  in  resistance  between  the 
two  probe  circuits.  Whereas,  there  are  resistance  differences  in  the  parts  of  the  circuit 
other  than  the  thin  wire  sensors,  there  is  a  bridge  balancing  circuit  that  automatically 
balances  slow  changes  in  the  resistance  (time  constant  ~  20s),  and  the  amplifier  filters  out 
changes  that  are  slower  than  the  expected  turbulence.  Therefore  it  is  the  difference  in 
resistance  of  the  two  fine  wire  sensors  that  is  actually  analyzed.  The  difference  in 
resistance  in  the  two  sensors  is  then: 


AR  =  Rl(Tl)-R2(T2 ) 

=  R0l+  aR01  (7j  -  T0)  -  R()1  - aR02 ( T2  -  T0) 


(12) 
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Further,  define: 


*0  -“(-^01  +^02) 

$  =  R<n  ~R02 

Substituting  the  definitions  into  Equation  (7)  and  rearranging,  we  obtain: 


AR  =  aR0(Tl-T2)  +  8 


(13) 


In  practice,  the  reference  temperature  should  not  be  too  different  from  the  measurement 
temperatures,  consistent  with  the  linear  approximation,  and  a  is  a  small  number,  so  the 
tenn  in  brackets  is  usually  not  much  larger  than  1 .  If  the  probe  pairs  are  chosen  to  have 
nearly  identical  resistances  such  that  8  «  R0  then  the  final  working  equation  is: 


AR  =  aR0AT  (14) 

where  AT  is  the  difference  in  the  temperature  of  the  two  probe  sensors,  which  is  the 
basis  for  the  structure  function  calculation  that  defines  the  turbulence  level. 

4.4  Calibration  Procedure 

The  Wheatstone  Bridge  of  the  thennosonde  is  driven  with  a  small  AC  current  alternating 
at  3kHz.  The  wires  operate  in  low  overheat  ratio  mode,  or  as  cold  wire  probes.  Because  of 
the  dynamic  bridge  rebalancing  and  the  need  to  reduce  noise,  the  thennosonde  circuit 
includes  a  passband  filter  to  eliminate  DC  and  high  frequency  noise  from  the  turbulence 
signal:  The  bandwidth  is  determined  by  a  4  pole  Butterworth  and  low  pass  filter  network 
whose  3dB  frequencies  are  1kHz  and  0.2Hz  respectively.  The  turbulence  signal  in  this 
passband  is  suitably  amplified  and  finally  passed  through  an  RMS  module  with  a  nominal 
time  constant  of  3.75s.  The  RMS  voltage  is  the  thennosonde  output  that  is  fed  into  the 
radiosonde  telemetry  system. 

As  a  result  of  Equation  9,  above,  a  linear  relation  is  expected  between  the  input  AT  and 
the  output  Vrms;  however,  this  is  not  assumed.  Since  the  board  filters  out  all  but  the  narrow 
turbulence  and,  a  special  calibration  instrument  was  designed  and  fabricated  that  switches 
known  resistors  in  and  out  of  the  Wheatstone  bridge,  in  place  of  the  wire  probes.  The  two 
probe  wires  are  brought  into  the  box  and  connected  to  two  nearly  identical  resistors  of 
27  Q  each.  One  of  those  resistors  is  part  of  a  parallel  circuit  that  is  designed  to  allow  one 
or  more  other  known  resistors  to  be  switched  in  and  out  of  the  circuit  at  20Hz.  In  other 
words,  one  probe  circuit  alternates  between  a  resistance  27  Q  and  a  lower  resistance  based 


8 


on  the  value  of  the  switched  in  resistor.  The  calibration  box  effectively  provides  a  square 
wave  input  of  resistance  difference  to  one  leg  of  the  bridge,  varying  between  0  and  a 
specified  resistance  value,  A R  .  The  bridge  balancing  circuit  balances  the  two  sides  so 
that  this  square  wave  becomes  symmetrically  positioned  about  0  volts. 

The  nominal  temperature  disturbances  are  assumed  to  be  either  sinusoidal  or  linear 
combinations  of  sinusoidal  disturbances.  Since  root  mean  square  of  a  square  wave  is 
twice  that  of  a  sine  wave,  the  resistance  is  halved  to  give  the  equivalent  of  a  sinusoidal 
disturbance.  Hence,  for  the  calibration  process,  the  temperature  difference,  input 
calibration  resistance  difference  relationship  is  detennined  by  the  formula: 


AT,  =  AR(:al 


'-cal 


2  a 


(15) 


Each  board  is  calibrated  by  measuring  the  voltage  of  the  output  of  the  RMS  chip  on  the 
board  for  a  progression  of  input  resistance  differences.  Note  that  while  a  is  a  property  of 
the  probe  wire  material,  R0  can  vary  since  it  is  also  a  function  of  probe  wire  geometry. 

Even  though  the  calibration  board  has  nominal  resistance  ( R0 )  of  27  Q ,  each  probe  pair 
will  have  an  R0  that  is  the  mean  of  the  reference  resistance  of  each  probe  pair,  and  will 
vary  about  27  fi  .  Therefore,  the  AT  •  R0  product  is  recorded,  which  is  ARcal  /  2 a  . 


4.5  Calibration  Example 

The  calibration  procedure  begins  by  attaching  the  leads  that  would  normally  go  to  the 
probes  to  the  calibration  box  inputs.  After  checking  to  insure  proper  input  current  to  the 
bridge,  a  DC  Voltmeter  is  attached  to  the  output  of  the  RMS  chip,  and  the  calibration  box 
is  set  to  the  maximum  A R  value  of  0.052234  Q  .  The  gain  of  the  amplifier  on  the 
thennosonde  board  is  set  to  produce  an  output  of  1.950V  out  of  the  RMS  chip.  Next  A R 
is  set  to  0,  and  the  board  is  allowed  to  settle  to  the  noise  floor,  typically  around  0.02V, 
which  is  recorded.  An  example  of  a  typical  calibration  sheet  is  shown  in  Table  1.  Column 

1  is  the  input  AR  ;  column  4  is  the  output  voltage  from  the  RMS  chip.  The  input  resistance 
is  set  for  each  of  the  remaining  values  of  AR  and  the  voltages  are  recorded.  Also  shown 
in  the  table  are  some  ancillary  calculations:  Column  2  shows  the  input  AR  divided  by 

2  a  ,  column  2  is  the  equivalent  AT  for  the  nominal  R0  of  27  Q  .  Finally,  the  last  column 
is  the  voltage  output  after  the  noise  correction,  Equation  5,  for  each  input  resistance  value. 
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DELTA  R 

Delta  R/2alpha 

delT 

Board  265 

265  COR 

0 

0.021 

0.000 

0.001042 

0.138933 

0.005069 

0.046 

0.041 

0.003137 

0.418239 

0.015259 

0.121 

0.119 

0.005215 

0.695330 

0.025368 

0.196 

0.195 

0.010445 

1.392715 

0.050810 

0.392 

0.391 

0.015676 

2.090128 

0.076254 

0.586 

0.586 

0.020894 

2.785918 

0.101639 

0.785 

0.785 

0.025670 

3.422729 

0.124872 

0.964 

0.964 

0.031088 

4.145021 

0.151223 

1.166 

1.166 

0.036600 

4.880032 

0.178038 

1.372 

1.372 

0.041801 

5.573507 

0.203338 

1.563 

1.563 

0.047018 

6.269031 

0.228713 

1.757 

1.757 

0.052234 

6.964583 

0.254089 

1.950 

1.950 

Table  1.  Calibration  sheet  for  thermosonde  board  256,  which  is  a  typical  example. 


4.6  End  to  End  Check  of  the  System 

The  above  calibration  procedure  insures  that  the  temperature  difference  sensed  by  the 
resistance  wires  is  properly  converted  to  volts  at  the  output  of  the  RMS  chip.  This  is  only 
the  first  of  the  several  operations  perfonned  by  the  thermosonde  system  before  data  is 
output  to  the  user.  The  entire  process  is  shown  in  Figure  4.  In  calibration,  the 
thennosonde  signal  processing  board  is  checked  for  volts.  When  fully  assembled,  that 
voltage  is  sent  to  the  radiosonde  interface  and  analog-to-digital  unit.  The  digital  unit 
converts  volts  to  counts  and  sends  the  signal  into  the  radiosonde  where  it  is  send  down  to 
the  ground  station  as  spare  channel  information  with  the  regular  meteorological 
information.  That  signal  is  sent  on  a  403MHz  carrier  to  the  receiver  of  the  Ground 
Station.  The  received  temperature  difference  data  is  processed  then  output  as  counts.  The 
counts  are  then  converted  back  to  voltage  by  post-processing  software,  and  then  further 
converted  to  temperature  difference.  These  conversions  from  volts  to  counts  to  radio 
signal  then  back  again  involve  two  conversion  factors  which  are  unique  to  each  interface 
unit  and  radiosonde.  An  end  to  end  check  of  the  system  was  perfonned  in  the  laboratory 
where  the  voltage  output  from  the  thennosonde  board  was  compared  to  the  voltage  output 
from  the  data  reduction  software.  The  results,  shown  in  Table  2,  are  quite  encouraging,  in 
that  quantities  are  either  identical,  or  differ  by  1  or  2  in  the  4th  significant  figure. 


Volts  Out 

1.154 

1.389 

1.187 

0.967 

0.787 

0.586 

0.197 

Counts  to  Volts 

1.153 

1.391 

1.186 

0.967 

0.787 

0.586 

0.197 

Table  2.  Comparison  of  voltage  from  board  to  voltage  from  the  data  reduction  program. 
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Figure  3.  The  probe  calibration  box  is  used  to  simulate  the  probes  as  described  in  Section 
4.5.  As  the  AR  is  switched  over  its  widest  range,  the  RMS  voltage  output  from  the  signal 
conditioning  electronics  is  compared  to  the  RMS  voltage  obtained  through  the  whole 
system  as  would  be  in  an  actual  flight. 


4.7  Universal  Calibration  Curve  Fits 

With  improved  quality  control  of  thermosonde  parts  and  assembly,  it  was 
observed  that  after  adjusting  the  gain  of  a  board  at  the  highest  calibration  setting,  the 
output  voltage  at  intennediate  settings  appeared  to  be  very  reproducible.  Following  the 
calibration  analysis  reported  in  this  section,  it  was  determined  that  a  single  (universal) 
calibration  curve  could  be  used  for  the  thennosondes.  The  single  universal  calibration 
curve  was  the  result  of  the  calibration  of  16  thermosonde  boards. 

As  described  in  section  4.2,  the  new  procedure  assumes  that  the  input  to  the  RMS  chip  is 
composed  of  signal  and  noise.  With  this  assumption,  the  output  of  the  RMS  chip  is 
always  corrected  by  Equation  5  to  obtain  the  signal.  For  historical  reasons  the  voltage 
signal  may  -  or  may  not  be  referred  to  as  Vrms.  The  calibration  curve  is  applied  to  the 
signal  voltage  after  correction  for  noise.  This  technique  has  the  advantage  of  insuring  that 
all  the  corrected  calibration  curves  run  through  the  origin.  The  data  from  the  calibration 
of  16  boards  are  shown  in  Figure  4.  The  equivalent  R0  AT  is  on  the  ordinate;  the 

corrected  (signal  input)  RMS  voltage  is  on  the  abscissa.  A  straight-line  calibration  curve 
fits  the  data  quite  well.  The  universal  calibration  formula  is: 
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(16) 


R0  AT  =  3.572  AVin 

As  shown  on  the  curve,  this  correlation  has  a  correlation  coefficient,  R,  of  1,  which  is  the 
result  of  rounding  up  the  value:  0.999969. 


Calibration  Curve 


RMS  Voltage 


Figure  4.  Calibration  points  from  16  thennosonde  boards  plotted  with  universal 
calibration  curve. 


4.8  Uncertainty  due  to  the  Universal  Calibration 

After  achieving  the  above  curve  fit,  uncertainty  estimates  were  made  by  comparing  the  fit 
to  each  of  the  16  calibration  points  for  each  Delta  T*Ro  level.  For  each  level,  the 
uncertainty,  as  defined  as  the  variance  of  the  differences  from  the  calibration  curve  and  the 
actual  points,  was  computed.  For  purposes  of  comparison,  a  nominal  value  of  Ro  equal  to 
27f2  is  used  to  depict  the  error  in  terms  of  temperature  uncertainty.  The  results  of  the 
calculations  are  shown  in  Table  1.  A  plot  of  the  difference  between  calibration  curve  and 
calibration  points  for  16  boards  is  shown  in  Figure  6.  The  maximum  errors  are  less  than 
10%  of  the  measured  value,  which  occurs  at  the  lower  temperature  value,  as  shown  in 
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Figure  7.  For  a  temperature  difference  of  0.0 15K,  the  maximum  error  is  around  2%  or 
less,  with  mean  errors  less  than  a  percent.  The  standard  deviation  of  the  fit  from  the 
calibration  points  was  used  as  the  uncertainty,  which  has  a  maximum  of  ~6%  at  the  lowest 
temperature,  and  decreases  1.5%  at  0.0 15K,  to  well  below  a  percent  as  temperature 
increases. 


Del  T*Ro 

Del  T  (K) 

Avg  Err 

%Err 

Var (K) 

St  Dev  of  Fit 

St  Dev  (%) 

Var  of  Fit 

0 

0 

0 

0.00% 

0 

0.00E+00 

0.00% 

0.00E+00 

0.138933 

0.005146 

-0.000284 

-5.52% 

3.92E-07 

3.17E-04 

6.16% 

1.01E-07 

0.418239 

0.01549 

-0.00011 

-0.71% 

6.59E-07 

1.93E-04 

1 .25% 

3.74E-08 

0.69533 

0.025753 

3.57E-05 

0.14% 

1.15E-06 

2.10E-04 

0.81% 

4.40E-08 

1.392715 

0.051582 

0.000197 

0.38% 

2.82E-06 

3.82E-04 

0.74% 

1.46E-07 

2.090128 

0.077412 

0.000139 

0.18% 

5.04E-06 

4.55E-04 

0.59% 

2.07E-07 

2.785918 

0.103182 

-1.92E-06 

0.00% 

6.04E-06 

4.73E-04 

0.46% 

2.24E-07 

3.422729 

0.126768 

-0.000165 

-0.13% 

8.27E-06 

5.79E-04 

0.46% 

3.35E-07 

4.145021 

0.153519 

-0.000244 

-0.16% 

9.9E-06 

6.56E-04 

0.43% 

4.30E-07 

4.880032 

0.180742 

-0.000345 

-0.19% 

5.99E-06 

5.91  E-04 

0.33% 

3.49E-07 

5.573507 

0.206426 

0.000192 

0.09% 

6.36E-06 

5.30E-04 

0.23% 

2.81  E-07 

6.269031 

0.232186 

0.000331 

0.14% 

9.56E-06 

6.93E-04 

0.30% 

4.80E-07 

6.964583 

0.257948 

-1.14E-05 

0.00% 

1.78E-10 

1.21E-05 

0.00% 

1.46E-10 

Table  3.  Comparison  of  the  universal  calibration  curve  to  the  calibration  values  from  the 


16  boards. 


As  a  result  of  the  new  universal  calibration  curve,  the  calibration  procedure  now  consists 
of  setting  the  highest  voltage  and  checking  intermediate  levels  to  insure  that  the  boards 
conform  to  past  standards  and  that  the  universal  curve  is  applicable. 
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Calibration  Curve  Error 


Delta  T  (K) 


Figure  5.  Difference  between  universal  calibration  curve  and  calibration  points  expressed 
as  Temperature  difference  assuming  a  nominal  resistance  value  of  210.. 


5.  INSTRUMENTATION  ERRORS 

The  thermosonde  instrument  is  composed  of  standard  meteorological  radiosonde 
components  and  the  probes  temperature  fluctuation  conditioning  circuit.  Instrumentation 
errors  in  optical  turbulence  measurements  can  be  attributed  to  errors  in  the  temperature 
and  pressure  sensors  of  the  radiosonde,  and  temperature  fluctuation  errors  in  the 
thermosonde  measurements. 

5.1  Pressure  and  Mean  Temperature 

The  pressure  sensors  of  current  radiosondes  have  an  accuracy  of  0.5  mbar,  which  is 
0.05%  of  sea  level  pressure  to  5%  at  maximum  altitude.  The  temperature  sensor  has  an 
accuracy  of  0.2K  at  sea  level,  which  degrades  to  0.4K  at  maximum  altitude.  This  is,  at 
worst  0.2%  of  the  ambient  temperature  throughout  the  range10.  Systematic  errors  are 
minimized  by  calibration  of  the  pressure  sensor  of  each  radiosonde  against  a  high 
precision  barometer  at  launch. 
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Error  as  percent  of  value 


Figure  6.  Error  between  universal  curve  and  calibration  points,  as  percent  of  the 
temperature  difference.  The  mean  error  at  any  input  temperature  is  shown  as  connected 
squares.  Finally,  the  standard  deviation  of  the  calibration  from  the  curve  fit  is  shown  as 
the  line  which  remains  in  the  positive  portion  of  the  graph. 


5.2  Temperature  Fluctuation 

The  thermosonde  measures  temperature  fluctuations  that  are  converted  to  optical 
turbulence  levels. 

5.2.1  Accuracy 

During  calibration,  voltages  are  read  within  lmV  for  a  range  of  voltages  from  the 
minimum,  which  is  the  noise  floor  of  about  lOmV  (equivalent  to  about  0.0005K)  up  to  a 
maximum  output  of  2.0V1 1 
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5.2.2  Precision 


With  lmV  accuracy,  the  precision  is  two  significant  figures  for  noise  floor  voltage 
readings  up  to  more  than  three  significant  figures  at  the  highest  output  voltages. 

5.2.3  Universal  Calibration  Curve 

As  shown  above,  the  universal  curve  causes  a  small  error,  which  has  an  uncertainty  of 
about  6%  at  0.005K,  to  well  below  a  percent  at  values  above  0.02K. 

5.2.4  Noise  Floor 

The  magnitude  of  the  noise  level  of  the  current  thermosonde  has  benefited  from  the 
aforementioned  improvements.  The  noise  level  has  decreased  from  an  original  level  “ 
equivalent  to  0.004K,  to  approximately  0.002K  in  the  1980’s13,  to  less  than  0.001K 
today. 

5.3  Effect  of  the  Above  Errors 

2  9 

Any  calculation  of  Cn  based  on  measured  Ct  is  subject  to  uncertainties  in  the  measured 
pressure,  temperature,  and  temperature  fluctuation.  Assuming  that  these  errors  are 
uncorrelated  with  each  other,  standard  uncertainty  analysis  of  the  Equation  (1-9)  results 
in  the  following  estimate  of  the  error  in  Cn2. 


where  cr,  is  the  total  error  and  uP ,  aT  ,and  <jst  are  the  uncertainties  in  pressure, 

temperature,  and  temperature  fluctuation,  the  square  root  of  Ct  .  The  measured 
temperature  differences  vary  from  the  noise  floor  to  within  a  tenth  of  a  degree.  The 
pressure  uncertainty  factor  can  become  important  at  the  highest  altitudes  for  high  values 
of  Ct2.  Temperature  uncertainty  is  never  significant.  For  the  lowest  values  of  Cj2,  the 
error  is  less  than  10%,  and  decreases  to  well  below  a  percent  as  Cfi  increases. 

For  a  modest  Cj2  value  of  0.0002  K2nT2  3,  the  uncertainty  in  Cn2  is  down  to  1.5%.  By  a 
stronger  Ct  of  0.01,  overall  Cn“  uncertainty  is  below  0.5%  at  low  altitudes,  and  starts  to 
noticeably  increase  at  20km  reaching  a  4%  uncertainty  at  30km  due  to  pressure 
uncertainties. 
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5.4  Other  Instrument  Related  Factors: 


5.4.1  Heat  Transfer  -  Altitude  Limits 

As  the  thermosonde  approaches  the  maximum  altitude  of  30km,  the  heat  transfer 
coefficient  has  decreased  to  the  point  that  the  cutoff  frequency  of  the  instrument  response 
approaches  sizes  that  can  impact  the  lm-structure  function.  Therefore,  altitudes  higher 
than  30km  are  not  recommended. 

5.4.2  RMS  Chip 

The  proper  computation  of  optical  turbulence  requires  an  accurate  running  root  mean 
square  of  the  temperature  differences.  In  fact,  commercially  available  RMS  chips 
perform  only  an  approximation  to  a  true  root  mean  square.  The  effective  time  constant  of 
the  calculation  is  not  constant.  When  the  temperature  differences  are  short  relative  to  the 
~ls  sampling  period,  which  is  associated  with  very  low  values  of  integrated  turbulence, 
there  can  be  an  error  of  a  factor  of  2  in  integrated  optical  turbulence.  Numerical 
simulations  with  more  typical  distributions  of  simulated  turbulence  show  errors  of  10  to 
20%  in  integrated  turbulence8. 

5.4.3  Telemetry  system 

The  output  from  the  RMS  chip  is  routed  to  an  interface  card  with  an  accuracy  of  less  than 
2mV,  which  is  well  below  the  other  errors  mentioned  in  the  report. 

5.4.4  Environmental  Factors 

Thermosondes  normally  ascend  at  approximately  5m/s  from  ground  level  to  30km 
attached  to  a  large  weather  balloon  by  a  1 10m  line.  There  are  some  inherent  sources  of 
uncertainty  that  must  be  considered:  the  wake  of  the  balloon  and  inherent  statistical 
issues  in  characterizing  turbulence,  with  its  random  variations  in  space  and  time. 

5.4.5  Balloon  Wake 

The  impact  of  the  balloon  wake  is  currently  under  investigation.  A  comparison  of  output 
from  a  normally  ascending  thennosonde  to  that  of  a  downward  facing  thermosonde  on  a 
descending  balloon  showed  that  wake  encounters  have  an  insignificant  effect  in  the 
troposphere,  but  can  be  significant  in  the  stratosphere14.  Encounters  with  a  wake  can 
occur  in  low  shear  environments15  when  measured  values  of  atmospheric  turbulence  are 
usually  low. 
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5.4.6  Stationarity  -  Homogeneity  Issues 

Since  the  earliest  thermosonde  measurements  (References  6,  7,  and  13),  there  has  been  a 
concern  about  the  ability  to  achieve  an  accurate  mean  value  of  the  temperature 
fluctuations  while  ascending  at  around  5m/s.  This  is  of  particular  concern  when  the 
layers  are  quite  thin,  since  the  instrument  is  only  sampling  the  RMS  chip  at  a  frequency 
of  about  1Hz.  The  turbulence  layers  that  cause  the  most  significant  decrease  in  optical 
performance  are  usually  thick  enough  that  many  samples  are  indicating  nearly  the  same 
value.  However,  accurate  measurement  of  thinner  layers  remains  an  unquantified  source 
of  uncertainty.  A  good  statistical  estimate  of  the  mean  requires  sampling  over  several 
correlation  lengths.  The  correlation  length  of  the  thennosonde  measurement  is  roughly 
equivalent  to  the  largest  scales  contributing  to  a  1  meter  structure  function  and  this  is  of 
the  order  of  10-20  meters.  Hence,  the  intrinsic  difficulty  in  obtaining  good  statistical 
estimation  of  Cn  on  an  ascending  balloon  in  layers  less  than  100  m  thick  is 
acknowledged. 

6.  Summary 

In  summary,  the  thennosonde  has  the  most  uncertainty  when  measuring  low  values  of 
turbulence.  Near  the  noise  floor,  errors  approach  the  size  of  the  measurement.  As  the 
turbulence  levels  increase,  uncertainties  from  all  sources  decrease  to  below  10%  as  the 
integrated  turbulence  levels  approach  that  of  Clear  1 .  The  strongest  levels  of  optical 
turbulence,  which  are  the  most  important,  reflect  the  region  of  highest  accuracy  of  the 
instrument. 
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